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Abstract —We consider the detection of a correlated random 
process immersed in noise in a wireless sensor network. Each 
node has an individual energy constraint and the communication 
with the processing central units are affected hy the path loss 
propagation effect. Guided by energy efficiency concerns, we 
consider the partition of the whole network into clusters, each 
one with a coordination node or cluster head. Thus, the nodes 
transmit their measnrements to the corresponding cluster heads, 
which after some processing, communicate a snmmary of the 
received information to the fusion center, which takes the final 
decision about the state of the nature. As the network has a 
fixed size, communication within smaller clusters will be less 
affected by the path loss effect, reducing energy consumption 
in the information exchange process between nodes and cluster 
heads. However, this limits the capability of the network of 
beneficially exploiting the spatial correlation of the process, 
specially when the spatial correlation coherence of the process 
is of the same scale as the clusters size. Therefore, a trade-off 
is established between the energy efficiency and the beneficial 
use of spatial correlation. The study of this trade-off is the main 
goal of this paper. We derive tight approximations of the false 
alarm and miss-detection error probabilities under the Neyman- 
Pearson framework for the above scenario. We also consider 
the application of these results to a particular network and 
correlation model obtaining closed form expressions. Finally, we 
validate the resnlts for more general network and correlation 
models through numerical simulations. 

Index Terms —Distributed detection; correlated measurements; 
energy and bandwidth constraints; wireless sensor networks; 
hierarchical clustering; spectrum sensing; cognitive radio. 

I. Introduction 

D istributed detection based on wireless sensor net¬ 
works (WSN) is a topic which has attracted great interest 
in recent years (see Q and references therein). A typical WSN 
has a large number of sensor nodes which are generally low- 
cost battery-powered devices with limited sensing, computing, 
and communication capabilities. Sensors acquire noisy mea¬ 
surements, perform simple data processing and propagate the 
information into the WSN to reach a decision about a physical 
phenomenon occurring in the coverage area. 

Network resources, such as energy and bandwidth, are 
scarce, expensive, and key variables when the design is 
focused on processing latency and detection performance. 
Clever data processing strategies are required to maximize 
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performance under resource constraints. Appropriate choices 
of the processing strategy could largely impact on the total 
cost or on the life cycle of the WSN, and as such, make its 
deployment on remote locations economically viable or not. 

A. Related work 

Distributed detection theory has been much studied in the 
past. Starting with the seminal work of Tenney and Sandell Q, 
several results have been derived on how each node process 
the available information and communicates with the fusion 
center (FC) where the final decision on the true state of nature 
is taken. Under this setup, digital transmission schemes, where 
appropriate quantization rules have to be designed have been 
considered in 0-0 (see also references therein). On the other 
hand, analog communication schemes were also studied in the 
past (see for example 

Instead of independent observations, often assumed in the 
previously mentioned works, we explore a more realistic 
scenario with spatial dependence measurements. Distributed 
detection problems with correlated observations can be consid¬ 
erably more challenging than their conditionally independent 
counterparts GD- Moreover, new theoretical frameworks have 
been proposed recently p3| , p4| . 

The design of distributed processing strategies exploiting 
spatial and/or temporal correlation among data is, in general, 
an open problem. It is well known that signal correlation can 
help to improve the detection performance, specially when 
low quality sensor measurements are available 0, Og, fig. 
Clearly, a clever use of the correlation in space and/or time 
requires cooperation among nodes. This problem has also 
been studied in previous works GD-Gg. Additional studies 
have explored the effects of correlation on the performance of 
distributed detection systems |2g, 1^. 

An important application that has gained a lot of attention 
in the recent years under this scenario is spectrum sensing 
for cognitive radio p2[ , p3) . Here, unlicensed, or secondary 
users (SU), want to detect the presence or absence of the 
licensed, or primary users (PU), in order to use the spec¬ 
trum when it is available. SUs could cooperate and build a 
WSN to achieve this task. When a PU uses the spectrum, 
it transmits a signal that may be sensed by the WSN. Due 
to the shadowing propagation effect, different close nodes 
sense correlated measurements. Thus, dependent observations 
must be considered in order to develop reliable spectrum 
sensing schemes. Moreover, the signal model and the scheme 
developed in this work can be directly applied to a cognitive 
radio scenario GD- 


2 


Typically, the nodes of a WSN are distributed on a vast area 
to sense the whole region of interest. Thus, the energy needed 
for reliable communication increases with the transmision 
distance due to the effect of path loss. In order to save 
transmission energy, it is convenient to reduce the region 
where sensors transmit high precision data. In this way, a 
hierarchical WSN is introduced; sensors are grouped in clus¬ 
ters with a designated cluster head (CH). Each CH processes 
the data independently of the other ones and communicates a 
summary of the information gathered by its nodes to the FC, 
where the final decision about the state of the nature is made. 
Previous works have already considered the idea of 

partitioning the region to be sensed into clusters. However, 
in p4) the authors considered conditionally independent and 
identically distributed measurements, and in Q, spatially 
coherence regions were assumed, but dependence among them 
is not taken into account. 


regime, where the energy dedicated to transmit the data 
limits the performance, and a correlation-limited regime, 
where the performance is limited by the poor use of 
the spatial correlation of the process. The ranges of 
validity of each of these regimes are given by different 
parameters of the WSN some of which are under control 
of the designer. The characterization of this problem 
gives crucial information for the deployment of a WSN. 

2) In some cases of interest the use of local correlation 
information is almost as good as the use of the global 
correlation statistics of the sensed processes. This is 
a very important issue from an implementation point 
of view, as the use of global correlation information 
requires a greater degree of complexity in the network. 

3) The use of local correlation information becomes more 
profitable when the sensor measurements are less reli¬ 
able (i.e., the measurement signal to noise ratio is low). 


B. Contributions 

In this paper, we develop a clustering WSN model with 
energy and bandwidth constraints. We assume that each cluster 
processes the information independently of other clusters. 
Moreover, we consider that each individual cluster only knows 
the local correlation structure of the process. This hypothesis 
simplifies the scheme in the sense that communication of high 
data rate information between CHs or nodes from different 
clusters is avoided. As these exchanges would require to over¬ 
come path loss effects, an energy saving in each node will be 
obtained. However, in a WSN with correlated measurements, 
this strategy limits the ability of the network of beneficially 
exploiting the correlation of the process. Therefore, a trade-off 
is established between the beneficial use of the underlying ran¬ 
dom process spatial correlation and the energy spent by each 
node in the network. This is a trade-off between the energy 
required to reliably communicate each node measurement and 
the ability to effectively use the cooperation scheme required 
for exploiting the spatial correlation information. 

To characterize this trade-off, we derive tight approxima¬ 
tions of the false alarm and miss-detection error probabilities 
for statistics built with sum of independent although not neces¬ 
sarily identically distributed random variables. The expressions 
involve the sum of the logarithmic moment generating function 
(LMGF) of each random variable and its first and second 
derivatives. For example, in this paper we consider Gaussian 
processes for which the statistic is a quadratic form of a 
Gaussian correlated random vector. In this case, the statistic 
can be expressed as a weighted sum of i.i.d. random variables 
and the error probabilities can be easily computed. 

We obtain closed forms expressions for a particular model 
of the sensed process. Also, we numerically compute the 
average performance for more general processes and when the 
sensors are randomly placed in space and show that similar 
conclusions can be drawn. These conclusions are summarized 
as follows; 

1) In a WSN with Gaussian measurements there are two 
clearly different operating regimes; an energy limited 


C. Notation and Organization 

Scalars are written in italic, vectors in boldface and matrices 
in capital letters. A„ and Bnm are square and rectangular 
matrices of sizes n x n and n x m, respectively. We define 
column vectors with components as a = where k 

belongs to the set of indexes A. Similarly, we denote matrices 
Anm = [(ik]k^A built through the horizontal concatenation 
of the column vectors {a^}, where G and the 

cardinality of A is |,A| = m. det(A„) and tr(A„) are 
the determinant and the trace of A„, respectively. (•)^ and 
(•)^ denote transpose and transpose conjugate. We do not 
distinguish between random variables and their realization 
values when the context is clear. p{x\'H) is the probability 
density function (p.d.f.) of x under H. Pi(') denote 

probability and the expectation respectively, both computed 
under hypothesis Hi- Given two statistics T and T', we denote 
that they have the same error probabilities writing T = T'. 
Given a function / 2 (-), its first and second derivatives are 
denoted as /i(-) and The notation f{n) = O {g{n)) 

means that \imn-^oo fin)/g{n) = c, where c is a non-zero 
constant. 

The paper is organized as follows. In Section]^ we present 
the WSN model and the binary hypothesis testing problem. 
In Section [I^ we describe the precoding strategies used for 
stationary and non-stationary random processes. In Section 


IV we derive a new global statistic and formulate the main 


problem of the paper. In Section]^ we compute tight approx¬ 
imations of the error probabilities to be used in Section |Vl] to 
analyse the performance of the proposed strategy. Finally, in 
Section [Vn] we elaborate on the main conclusions. Technical 
proofs are provided in the appendices. 


H. WSN Model 

We analyse a sensor network with three hierarchies; nodes, 
clusters with their corresponding CH, and the FC. The nodes 
are grouped in clusters. Each CH gathers information from the 
nodes in its corresponding cluster, performs some processing 
independently of other clusters and transmit a summary to the 
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FC, where the final decision is made. A sketch of the WSN 
is shown in Fig. [1] 



Fig. 1. A hierarchical WSN, where the information is propagated from the 
nodes to the CHs, and from them to the FC, which decides about the state of 
the nature. 


The WSN has m clusters and a total amount of n nodes 
which are arbitrary distributed in an area uAq, which for 
simplicity it will be assumed to be a square and Aq a constant 
value (with units). Each cluster has li nodes, i = 1,..., m, 
such that n = Let Ci be the set of sensors of the i- 

th cluster, i = 1,... ,to. {Ci}™ is a partition of the set of 
sensors [1 : n], i.e., Ci OCj = % ^ j, i = 1,..., m, and 

Ui^iCi = [1 : n]. In the following, we provide a description 
on the main elements of the general WSN model we will work 
with. 

The fc-th sensor of the WSN takes a measurement described 
under each hypothesis by 

f ’Hi : Vk = Sk+Vk, -j. 

\ T-l-o ■■ Vk =Vk, k=l,...,n, 

where Sk is a correlated zero-mean circularly-symmetric com¬ 
plex Gaussian process with variance and autocorrelation 
function (ACF) E[sjSfe] = Rs{j,k) with j,k S [1 : n]. 
In general, jsfe} is a non-stationary spatial process. We also 
assume that Vk is a zero-mean circularly-symmetric complex 
white Gaussian noise independent of Sk with variance a^. 
Thus, Uk is Gaussian distributed either under "Ho or "Hi. 
The above model can be easily generalized to the case in 
which each sensor takes several measurements at different time 
instants, and the processes have temporal correlation. All the 
mathematical developments and conclusions in this paper can 
be extended to that case. However, in order to keep things 
simple we consider only the case in which each sensor takes 
only one measurement. This is in line with several works on 
distributed detection on WSN ii-jTg. 

Define the vector = \yk]k<^Ci the set of measurements 
taken by the nodes in the z-th cluster. Define Si = [sk]k^Ci 
and Vi = [vk]k&Ci accordingly. The covariance matrices 
of Si and Vi are, respectively, Eg and where the 

elements of are Rs{j,k), with j,k G Ci, and R. is 

the identity matrix of dimension R. Therefore, the covariance 
matrix of under Hq is = cr^Ri, and under Hi is 

El,;, = Eg,;, +alR. Let y = [yiRe[i-.m] be the measurements 
taken by all nodes in the network and define s = 

V = [vi]i^[i-,rn\- The covariance matrices of s and v are, 
respectively, Eg,„ and where the elements of Es,„ are 
{T,s,n)j,k = Rs{j,k). The covariance matrix of y under Hq 


and Hi, are, respectively, 

LjO,n — ^l,n — Eg,n -f (Tyin- (2) 

Notice that Eo,n and Ei,„ can be thought as m x m block 
matrices, where the diagonal blocks are Eq,;, and Ei,;,, i = 


III. Cluster Precoding Strategies 

In this work, we consider that the communication between 
the nodes and their corresponding CH is through either a 
parallel-access channel (PAC) or multiple channel use^ of 
a multiple-access channel (MAC). In the case of the PAC, the 
amount of channel uses coincides with the number of available 
orthogonal channels in each cluster. In both cases, we assume 
that the nodes communicate with its corresponding CH without 
any interference from other clusters. Consider now that the i- 
th cluster has /' channel uses available for transmission, with 
I'i ^ k- Each channel use is associated to either an orthogonal 
time slot or a frequency band. Then, the processing strategy 
may use l[ time slots, Z' frequency bands or a combination of 
them such that the product of time slots and frequency bands 
is An illustration of the distributed scheme is shown in Eig. 

|2(i)] 



(a) 



(b) 


Fig. 2. A distributed scheme for gathering information in a WSN. (a) The 
nodes use a PAC or a MAC channel to communicate with their corresponding 
CH. (b) The FC collects the processed information of each CH and takes a 
decision. 


A. The PAC precoding strategy 

If a PAC is considered, the received signal at the ?-th CH 
is 

Zk = akCkVk+Wk, kGUi, (3) 

^We say that a channel use takes place when a symbol is transmitted. 
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where Wk is the zero-mean communication white noise with 
variance cr^ and circularly-symmetric complex Gaussian dis¬ 
tribution independent of the measurements yk, Ofc = ||£Cfc — 
xcHi\\~^ is the factor associated with the path loss effect 
between the position of the node G and the position 
of the corresponding cluster head aJcH; G e is the path 
loss exponent, is the precoding coefficient used by the A:-th 
sensor, and Ui C Ci is the set of channel uses indexes for the 
i-th cluster with cardinality \Ui \ = Z'. Under this strategy, each 
sensor has a single channel use to communicate its measure¬ 
ment. If only < li channel uses are available, li —1[ sensors 
in the cluster stay silent. Let Ci = d\2Lg{{ck)keUi) be the 
precoding matrix of the i-th cluster, Ai = dmg{{ak)keUi)-’ 
define the column vectors Zi = [zk]k^Ui, 
with covariance matrix The received signal vector in 

each cluster can be expressed as 

2 :^ = Ai C.Vi + i = 1,. .., m. (4) 

with covariance matrix under T-Lj, j = 0,1 given by, 

S,,,, = Cf . (5) 

The average energy consumed by the fc-th node in the 
corresponding i-th cluster during one of the scheduled 
transmissions is; 

= E (|cfe 2 /fcp) = +Pi(^l), (6) 

where pi is the a priori probability of the state of nature Hi- 
When Pi is unknown, a natural upper-bound for Ek is obtained 
by taking pi = 1. In this work, we will consider the following 
individual sensor energy constraint; 

Ek<Eyk€[l:n\. (7) 

For the particular case when the sensed process is uncorre¬ 
lated, and considering an individual sensor energy constraint, 
it can be easily proved that the optimal strategy consists in 
using all the available energy at each sensor. In this paper, 
we will consider this strategy even for correlated processes. 
The precoding strategy for a PAC channel named amplify and 
forward strategy (AFS-PAC) is defined next. 

Definition 1 (AFS-PAC). Considering the i-th cluster, (j^ and 
the precoding matrix for the AFS-PAC strategy is 

Ci = i = ( 8 ) 

B. The MAC precoding strategies 

In the case that a MAC is used, each sensor communicates 
with its corresponding CH using different gains for each chan¬ 
nel use. The signal collected at the i-th CH is a noisy version 
of the coherent superposition of the symbols transmitted by 
the li sensors through the corresponding MAC; 

Zk' = akCkk'Vk + Wk', k' e Ui, (9) 

keCi 

where Ok, Wk' and Ui are the same quantities defined above. 
Thus, the fc-th sensor transmits in the fc'-th channel use 


its measurement yk scaled by the precoding coefficient Ckk' 
through the MAC channel. Stacking the measurements from 
each cluster, the same vector signal model 0 is obtained, 
although for the MAC channel, the precoding matrix defined 
by Ci = \ckk']k^Ci,k'eUi is not, in general, a diagonal matrix. 
Considering (j^, the average energy consumed by the fc-th 
node in the corresponding i-th cluster during the Z' scheduled 
MAC channel uses is; 


FZfc = E 



^kk' Uk I 




k'eUi 


( 10 ) 


If {sfe} is spatial-stationary, two schemes presented in 114| 


may be used to optimally allocate the energy and bandwidth 
in successive uses of the MAC channel. Both schemes, called 
principal components strategy in a multiple-access channel 
(PCS-MAC) and principal frequency strategy in a multiple- 
access channel (PFS-MAC), are asymptotically the best strate¬ 
gies among the orthogonal schemes considering energy and 
bandwidth constraints and the error exponents as metrics of 
performance. Although error exponents are sometimes useful 
to compare and define strategies of large networks, they do 
not provide an exact characterization of the error probabilities, 
which are the typically required metrics of performance for 
WSNs. Therefore, in this work, we mainly use the error 
probabilities as metrics of performance. Because the optimum 
strategy for a MAC, in terms of error probabilities, remains 
an open problem, we will use the PFS-MAC strategy (which 
is optimal in terms of error exponents) and show that good 
results are obtained. 

On the other hand, if {s^} is not spatial-stationary, in 
general, the covariance matrix of the process does not have 
a regular structure (i.e. Toeplitz form). However, in this case 
the PCS-MAC scheme presented in p4l Def. 2] may be used 
to optimally allocate the transmission energy and bandwidth. 
This scheme requires a built infrastructure that might not 
available be available on any WSN. For that reason, in the 
spatial non-stationary case, we will consider the AFS-PAC 
strategy. 

Next, we describe the PFS-MAC strategy for stationary 
processes. During each channel use, each CH receives a noisy 
version of a given frequency bin of the discrete Fourier trans¬ 
form (DFT) of the corresponding measurement vector y^. The 
idea behind this strategy is to transmit to the corresponding 
CH the most distinguishing spectral components of the power 
spectral density (PSD) of the sensed process (j){i') (achieved 
when I grows unbounded). In particular, we will use for each 
cluster a simple almost-optimal solution that employs the PFS- 
MAC scheme with a constant energy profile (EP) on a given 
set of transmission frequencies. The strategy, called ON/OFF- 
EP (see |[T4]), is defined as follows; 

Definition 2 (PES-MAC with ON/OEE-EP). Consider the i-th 
cluster. Eor a given k, denote by {ji,j 2 , ■ ■ ■ ,jii) ^ permutation 
of {1,2,..., li\ such that f f > ■ ■ ■ > f ■ 

Let T be the number of channel uses or bandwidth constraint. 
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Let /3i = ^ be the fraction of channel uses (or bandwidth) 
used at cluster i. Then, the set of channel uses is dehned by 
Ui = {ji,..., j;'} and the precoding matrix of the PFS-MAC 
strategy with ON/OFF-EP is 


Ci = 7 iA- I = 1,... ,m, 


( 11 ) 


where Fi 


uv. 


= G C 


liXl', 


is a sub-matrix of 


the normalized DFT matrix of size li x li, i.e., <^j,/ = 
ifik', f 2 k',- --Jiik'V with fkk' = ^ exp 
k = 1 ,..., (i; the diagonal matrix A~^ performs a channel 
inversion to compensate for the different gains introduced by 
the path loss. Using ( 111 in ( [T 0 ] l, and considering the individual 
sensor energy constraint (j7)7we obtain that 7 ^ is limited by 
the longest distance between the nodes and the CH of the i-th 
cluster. Then, 


li = 


E 




, z = l,...,m, 


( 12 ) 


where dn,ax,z = 

Using the above scheme, the covariance matrices of Zi 
under Hi and Hq are, respectively. 


— li^hvJ^s-hFlil'. + iliCkl + 0'w)lv. 


(13) 

(14) 


IV. Detection Statistics 


A. Full Correlation Strategy 

Suppose that the FC has direct access to all the measure¬ 
ments available in the CHs 2 : = [zj,..., through ideal 
noiseless channels. The covariance matrix of 2 : under Hj, 
3 = 0,1, is 


“ ^nn'^j,nBn 


cr^ I 


(15) 


where = diag(AiC'i,..., AmUm). Consider the 

Neyman-Pearson problem p5] for a hxed false alarm prob¬ 
ability level a, where a false alarm event occurs when Hi is 
declared but Ho is true. Let the logarithmic likelihood ratio 
(LLR) p5) of the full correlation strategy (FCS) be; 


7 fcs( 2 :) = log 


p(2;|Ho) 




(16) 


where we discarded the constant term log det(So^„'S|j"^,) 
without affecting the performance of the statistic. Now, under 
the Neyman-Pearson setting, the optimal decision rule chooses 
Hi if Tfcs( 2 :) > T„, and Ho otherwise, where the threshold 
of the test r„ depends on a. 

In a distributed setting, with noisy links from the CHs to the 
FC, this statistic provides a lower bound for the false alarm 
and miss-detection error probabilities. It is worth to mention 
that these lower bounds are not necessarily tight due to the fact 
that they do not contemplate the degradation effect introduced 
by the communication channel between the CHs and the FC. 


B. Local Correlation Strategy 

When the process {sfc} is spatially correlated, the covari¬ 
ance matrix 5i in is neither diagonal nor block- 
diagonal. Hence, cannot be expressed as the sum of 
the local LLRs from each cluster, and processing the data 
independently does not lead to the global LLR. However, it is 
possible to implement simple distributed detection schemes to 
build an appropriate statistic at the FC using the compressed 
data transmitted by the CHs. 

The FCS statistic can be decomposed in such a way that 
each term is a function of one or at most two sets of 
measurements from each cluster, i.e., 

m 

Tfcs{z) — ^ 2 :^ ~ “ i . J ') 

m mm 

i=l i=l j = 


where 220 .;' ^ 1 are the (i, i)-th block of and 

respectively. Mu is the difference between the (i,i)-th block 


of 




and 




; and M^j is the (i,j)-th block of 


"l,n' ■ 


The first term in (17 1 captures most of the correlation between 


sensor measurements of the same cluster. We refer to this as 
the intra-cluster correlation term. The second term captures 
part of the intra-cluster correlation not included by the hrst 
term. The third term considers the correlation between sensor 
measurements of different clusters and we refer to this as the 
inter-cluster correlation term. 

Considering that the correlation between measurements is 
a local effect, it is expected that the inter-cluster correlation 
term does not affect too much the whole statistic if each cluster 
has an appropriate size respect to the coherence distance of 
the proces^ Keeping this effect in mind and considering 
also that gathering measurements from distant nodes is energy 
expensive, it would be reasonable to use a scheme that 
takes into account only the intra-cluster correlation. Then, 
the FC could processed the data sent by the CHs without 
considering the inter-cluster correlation. The strategy is named 
local correlation strategy (LCS). 

Let the local LLR at the i-th cluster be: 


TLCS,i{^^) = log 


p{z^\nl) 

p{z^\Ho) 


= Z 


H (-^-1 

\ 0,1 i 


(18) 


Each CH is able to exploit the intra-cluster correlation by 
computing the local LLR Ti^cs,i{zi) while keeping the com¬ 
munication energy cost relatively low. Moreover, no additional 
overhead is necessary to communicate and/or estimate the 
ACE between measurements of different clusters. Only the 
local intra-cluster ACE need to be known or estimated at each 
CH. In the next step of the scheme, each CH transmits its 
corresponding Tix:s^i{zi) to the EC. This communication can 
be done using a PAC or a MAC. When a PAC is used, since the 
inter-cluster correlation is discarded, the PC builds the global 


^Measurements taken in sensors separated by a distance greater than the 
coherence distance of the process are weakly correlated. 
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statistic adding the local statistics provided by the CHs, as 
suggested by the hrst term of ( [T7] i: 

m 

TlCs(^) = E (^ol - Sri;) (19) 

i=l 

On the other hand, if a MAC is used, the same statistic can be 
obtained. Notice, however, that the sum of the local statistics 
is naturally performed by sending them synchronously through 
the MAC channel. The channel communication between the 
CHs and the FC is assumed to be noiseless. We can justify 
this saying that usually there are much less CHs than nodes 
in the network and they typically have some infrastructure 
that the nodes do not have. For example, the CHs could 
have high directive antennas pointed to the FC, better energy 
budgets than the nodes, and, for the strategies considered 
in this paper, the CHs transmit scalar numbers which could 
be communicated using a reliable low rate coded digital 
modulation. For later developments, it is convenient to express 
the statistic as follows 

Tlcs{z) = z^T^,z. (20) 


D. Clustering partition problem 

Now, we are ready to formulate the main problem of this 
paper. We hrst summarize the whole strategy used. The n 
sensors of the WSN take measurements and transmit them to 
their corresponding CHs using the PFS-MAC strategy or the 
AFS-PAC strategy. Each CH has access to Q and allows the 
FC to build a statistic Tj using one of the previous strategies: 
j G {FCS, LCS, ED}. The EC makes a decision about the 
state of the nature comparing Tj against threshold tj. The 
false alarm and the miss-detection probability are, respectively, 
= Po(T, > T,) and = P,(T, < r,). 

The way each sensor node is assigned to each cluster is an 
important issue affecting the performance of the network. One 
may wonder which is the best allocation possible and which is 
its level of performance. To that end we can fix the probability 
of false alarm to a level a € (0,1) and consider an energy 
constraint given by Q. Let Ck be the label of the A:-th sensor 
that identihes to which cluster it belongs, i.e., = i if k 

belongs to the i-th cluster. We can then formulate a problem 
in order to obtain the optimal clustering partition that could 
be cast as follows: 


where r„/ = 


diag (^E 


'"-1 _ "-1 7^-1 _ 7 T -1 ^ k a 

"0,i( ■ ■ ■ ’^0,1'^ J 

n' X n' block diagonal matrix, with n' = Notice 

that, in general, r„/ ^ ^on' ~ Hence, the first term of 

( [20| i is not exactly the LLR of 2 ; although it could be a good 
approximation if the size of the clusters are greater than the 
coherence distance of the correlated process. In the EC, the 
statistic Tlcs(- 2) is compared against a threshold r„ in order 
to decide about the true state of nature. 

Now, we can elaborate an analytical justihcation about why 
we ignore the inter-cluster correlation in the LCS strategy 
using (17 I. The difference between the FCS and LCS statistics 
is the matrix used to produce the quadratic form: while Sj} 
is the same in both statistics because it is a diagonal matrix. 




“1 n' 7 ^ ’ ■ ■ ■ ’ )■ However, it can be shown, by 

using the Schur complement, that the Taylor approximation of 
order 0 of the l^ x l[ block matrices of is exactly the block 
diagonal matrix diag(5j";} ,..., Sj"}, ). Clearly, this would be 
a good approximation when the inter-cluster correlation is 
weak. 


C. Energy detector 

We also consider the statistic usually called energy detector, 
which does not contemplate the correlation of the sensed 
process. Therefore, it will result useful to compare it with 
the previous statistics to show how much one can loose if the 
correlation is not taken into account. Once each CH has the 
vector of measurements Zi in (j^, it builds T^'D^i{zi) = WziW"^ 
and sends it through a PAC or a MAC to the FC. The energy 
detector available in the FC is dehned by 

Ted(^) = ll^f. (21) 

Notice that this statistic is optimal only when the process {s^} 
is white noise process. 


inf P" 

(ei....,e„)G[l:m]" 

s.t. Pf" < a, Ek < P, Vfc G [1 : n]. 


( 22 ) 


where the threshold of the test tj depends on a. This is 
a challenging non-linear integer programing problem with 
constraints. In fact, it is a NP-hard problem and its solution is 
out of the scope of this paper. As our objective is to understand 
the trade-off between the benehcial use of spatial correlation 
against energy consumption in order to overcome the path-loss 
effect, we will consider a simpler clustering partition problem. 
We will assume that all the clusters have the same shape, in 
particular, we will assume that each of the m clusters are 
square regions of area liAg, where li is the number of sensors 
in the Pth cluster, and nAg = ^*^0 is total area in which 

the network is deployed. This simple model captures the basic 
fact that when the number of clusters is large, the size of them 
will be small, and the path loss effect in the local transmissions 
will be less harmful. However, as the clusters are smaller, less 
sensor nodes are present within them. In addition, some other 
nodes in their local neighborhoods, which can have strongly 
correlated measurements, could be assigned to other adjacent 
clusters. As the clusters do not cooperate and they make their 
processing independently, this clearly reduces the benefit of 
using the spatial correlation in each cluster. This basic model 
will be used in the following section with minor variations. For 
example, in Section fVI- A we will assume that each cluster has 
exactly I nodes arbitrary distributed in space. On the other 
hand, in Section |VI-B[ we will assume that the nodes are 
spatially distributed as a homogeneous Poisson Point Process 


(PPP) (261 where the average number of sensor per cluster is 


I, although for each realization of the spatial random process 
each cluster could have different number of nodes. In both 
cases, we variate the area of each cluster and, consequently, 
the number of sensors per cluster I (or its average) but keeping 
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constant the total amount of nodes n in the sensor network (or 
its average). 


V. Computation of the error probabilities 

In this section, we present a result used to compute tight 
approximations of the miss-detection probability P” and the 
false alarm probability using the LMGFs of the statistic 
under Hq and Hi. 

Theorem 1. Let {yk}k=i be mutually independent ran¬ 
dom variables with probability density function (PDF) pk 
and LMGF pfc(s) = logE Assume that E(?/^) and 

E(|j/fc — E(?/fc)p) exist and are finite Vfc S [1 : n]. Let 
Tn = yi + ■ • ■ + Vn with LMGF y,T„ — logE and 

let Tn S K. Then, if r„ > E(T„), 


P(T’« > Tn) = 

where sq > 0 is the solution to = Iit^{so). On the other 
hand, if < E(T„), 

P(T„ < Tn) = 

where si < 0 is the solution to r„ = /tT„(si). 

Proof: See App. |I-A| ■ 

Corollary 1 (Gaussian quadratic form). Consider the statistic 
given by the quadratic form 


P f — P / z 

±n' — n' ^ 5 


(25) 


where the distribution of jz under Hi is CA/'(0, 5^ „/), i = 0,1. 
Assume that and are positive definite matrices and 
that Tn' G (Eo(T„/),Ei(T„/)) is the threshold of the test. The 
false alarm and the miss error probabilities are, respectively. 


' .^27i-s2tr{[(Ho 


X e 


pn _ 




o 


X e 



1 X 

\VP)) 

I'—SO^Q^r 


(^)) 

1 X 

,'-Sl=l,r 



where sq > 0 and si < 0 are the solution to r„/ = 
tr{[(5o.n'T’n')“^ - soIn']~^} and r„- = ti-{[(Si_„/P„/)“i - 
respectively, and make /„/ — Si'E.i^n'Pn' positive 
definite for i = 0,1. 

Proof: See App. |I-B| ■ 


VI. Pereormance Analysis 

In Section |VI-A we consider a case of study where the 
process {s/c} is assumed to be equicorrelated, a particular 
case of a stationary process. It properly models the spatial 
correlation between the nodes and allows to obtain closed 


forms expressions of the error probabilities. In Section |VI-B| 
we will consider that the nodes are spatially distributed as a 
PPP We will see that similar conclusions to those obtained 
for the equicorrelated process can be drawn for non-stationary 
process. Before continuing, let SNRm = ^ and SNRc = ^ 
be the signal to noise ratio of the measurements and of the 
communication channel, respectively. 


A. A case of study: the equicorrelated process 

Consider a stationary zero-mean circular-symmetric Gaus¬ 
sian equicorrelated process with ACF Rs{k) = — p)Sko + 

alp, where 5ko is the delta of Kronecker. The PSD of the 
equicorrelated process is (t>{v) = a'l{l — p) + a1p5{v), where 
5{v) is the Dirac delta distribution and p is a scalar parameter 
that controls the correlation, 0 < p < 1. The covariance 
matrix for the measurements of each cluster is given by 
Ss,; = ~ P)^i + Tgpli where 1; is a Z x / matrix of 

ones. 

As the equicorrelated process is stationary, both PFS-MAC 
or AFS-PAC strategies can be used. However, in this section 
we will focus on the first one. We assume that all the clusters 
are identical and have the same parameters in Def|^ li = I, 

= I', Pi = P, di = (imax(0 ^nd 7 j = 7 o VL Considering that 
the CH could be in any position inside the coverage area IAq 
of its corresponding cluster, and that each cluster is a square 
region, the maximum possible distance between a node and 

its corresponding CH is dniax(Z) = \/2lAo. 

2 2 

Let SNRpFs-MAc(Z) = ..yp 2 ° 2+^2 be the effective signal-to- 
noise ratio of the measurements available in each CH. If we 
replace ([T^ in the previous definition, we have 


SNRpFS-MAc(0 


SNRmSNR c 

SNRc + (l-f SNRm)(2ZAo)5’ 


(26) 


where we have emphasized that the effective signal-to-noise 
ratio is a function of 1 . 

In order to evaluate the error probabilities of the LCS, FCS, 
and ED schemes, we need to compute the LMGFs of each 
statistic under Hi and their first and second derivatives to 
then apply Th. [T] For this particular process, it is possible 
to obtain closed forms solutions for the error probabilities for 
any threshold t„. However, these expressions are lengthy and 
are not shown here. The symbolic expressions are shown in 
the MATLAB script SymExprForPfaPm.m available with this 
paper. 

Therefore, we make focus on the case when the threshold 
for each strategy is set such that the miss error exponent is 
zero, where simple expressions are obtained. We summarize 
the results in the following lemma. 

Lemma 1 (Probabilities of false alarm for the equicorrelated 
process using the PFS-MAC strategy). Consider the following 
thresholds such that the miss error exponents are zero for the 
three strategies: rpcs = Ei(rLcs(2:)), tfcs = Ei(Tfcs( 2:)) and 
Ted = Ei(Ted( 2:)). Then, the false alarm probabilities of the 


















LCS, FCS and ED strategies using PFS-MAC are, respectively, 


P 

-ff 




fa,LCS \^ri/27rn[(i+/3-2)p2+2(l-/3)p+/3] 

g-n(r(p+^(i-p))-iiog(i+r(i+(;-i)p))-(/3-i)iog(i+r(i-p))) 


(27) 


o 




fa,FCS \^r^27rn[(n+/3-2)p2+2(l-/3)r+/3] 
„-n(r(p+/3(l-p))-ilog(l+r(l + (n-l)p))-(/3-i)log(l+r(l-p))) 


fa, ED 


( r ^ ( V" ) ) 


(28) 

„-n/3(r-iog(i+r))^ ^ 2 g-^ 


where F = SNRpfs-mac(0- 

Proof: See App. ■ 

It is easy to show that ^cs decreases monotonically 
with n, P and p. If SNRpfs-mac(0 were independent of 
I, Pf" LQg would also decrease monotonically with I and 
SNRpFs-MAc(0- However, as it is seen in (26 1 , SNRpfs-mac(0 


decreases with I, and therefore, a trade-off between these two 
parameters is established; for a given spatial density of sensors, 
to increase I, the number of nodes on each cluster, permits to 
successfully exploit the correlation, although this produces a 
worse SNRpFs-MAc(0 each cluster head. If the cluster size 
is reduced, SNRpfs-mac(0 improves but the correlation is not 
fully exploited. 

In Fig. the false alarm probability of the three schemes 
LCS, FCS and ED given by ( p7] l, ( |28| l and (291 respectively, 
are plotted against I, for the parameters shown in the caption. 
We observe that for the LCS scheme, to increase the size of 



Fig. 3. Trade-off between exploiting correlation and improving the measu¬ 
rement SNR. Pai'ameters: SNRc = 12 dB, SNRm = —6 dB, p = 0.9, e = 2, 
= 1 and n = 100. 

the cluster allows to exploit the correlation of the process and 
decrease Pfapcs- Beyond a certain cluster size, say I = 10, the 
benefit of exploiting the correlation is overcome by the lost 
of SNRpFs-MAc(0’ therefore, causing a worse performance. 
In the case of ED, increasing the cluster size reduces the 


effective SNRpfs-mac(0 without exploiting the correlation, 
due to the nature of this scheme. On the other hand, as the 
ECS scheme has access to the whole vector of measurements 
z, and computes the optimum LLR detector, it is only limited 
by SNRpFs-MAc((). which decreases when I increases. Eor a 
fixed cluster size, the ECS performance allows to quantify how 
much one can gain potentially by exploiting the correlation of 
the sensed process. Eor example, if I — 10, one could expect 
to decrease P^lcs between 3 and 4 orders of magnitude, 
although this is not entirely true given that the ECS scheme 
permits the EC to have access to the full vector z neither being 
penalized in bandwidth nor in transmitted energy, i.e., it could 
be a loose lower bound for some Z’s. Additionally, in the case 
of the PES-MAC scheme, the energy received in the CHs is 
limited by the furthest node, which could impose severe energy 
limitations. Notice that by definition of the LCS strategy, it 
coincides with the ED strategy when I = 1 (no correlation is 
considered in the statistic) and with ECS strategy when I = n 
(the network has a unique cluster). 


In Eig. 1^ we plot the false alarm probability of the three 
statistics in Lem. [T] against SNRc with I as a parameter. In Eig. 
|4(a)| we plot only the LCS strategy for the sake of clarity. Two 
regimes can be observed there; i) an energy-limited regime 
for which SNRc limits the performance and, ii) a correlation- 
limited regime, where the performance is limited due to fact 
that the inter-cluster correlation is not considered. When the 
cluster size is small, e.g. low number of sensors per cluster I, 
the sensors are near to the CH and the communication between 
them is reliable. We see in Eig. |4(a) that a WSN with small 
clusters performs better when the sensor energy budget is tight 
or the path loss effect is severe. However, when the sensor 
energy budget and the path loss effect are not an issue, large 
clusters are preferred because they can exploit successfully 
the intra-cluster correlation, obtaining large correlation gains 
even for weakly correlated processes (p = 0.2), as indicated 
in the figure. In Eig. 4(b) we show the behavior of the ECS 
and the ED strategies. We see that both strategies have lesser 
false alarm probabilities for small cluster sizes for the entire 
SNRc range. In other words, their performance is better due 
to an energy gain of the measurements in the CHs. Neither 
ECS nor ED modify how they consider the correlation of the 
process when the cluster size changes; the ED strategy always 
obviates the correlation while the ECS always considers the 
full correlation of the process (intra and inter-cluster). On the 
other hand, the proposed LCS makes the compromise between 
reliable measurement transmission to the CHs and exploiting 
the correlation of the process. Eor a given I, the ECS and 
the ED strategies provide a lower bound and an upper bound, 
respectively, of the error probability of the LCS strategy. 


Similarly to Lem. [T] we use Th. [T] to compute the miss- 
detection probability for the three statistics LCS, ECS and ED, 
when the threshold for each test r„ is such that the false alarm 
error exponent is zero, i.e., is the mean of each statistic 
under Hi. The expressions are lengthy and are not shown 
here. We have plotted them in Eig. against SNRc with I as 
a parameter. It is observed that, similarly to the probability of 
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Fig. 4. Probability of false alarm versus communication signal to noise ratio 
with parameters SNRm = — 3 dB, p = 0.2, e = 2, /3 = 1 and n = 100. 


false alarm, two regimes are observed for the miss-detection 
probability: a energy-limited regime and an correlation-limited 
regime. For a low SNRc, the miss-detection probability is 
limited by SNRpfs-mac(0- Therefore, a lower is obtained 
when the cluster size is small (low 1). For a high SNRc, 
the performance is limited by how much the correlation is 
exploited. As SNRpfs-mac(0 an issue, the LCS can 

increase the cluster size in order to exploit the intra-cluster 
correlation and decrease P^. 

To conclude this section, we obtain the optimum size Zopt 
which minimizes the false-alarm error probability of the LCS 
strategy given by (27 1 among all cluster sizes Z G [1 : n] such 
that n = ml. This is done numerically given that there is not a 
closed-form expression for Zopt- In Fig. 6(a) we plot Zopt versus 
SNRc with SNRm as a parameter. We see that the optimum 
cluster size increases with SNRc for each SNRm. We also 
observe that for a fixed SNRc, lopt is greater for lower SNRm. 
That is, as SNRm decreases, it is better to exploit the intra¬ 



Fig. 5. Pm versus SNRc- SNRm = 0 dB, n = 100 and p = 0.9. 


cluster correlation increasing I than to increase SNRpfs-mac(0- 
This is an important guideline for the case in which the WSN 
is composed by cheap nodes with limited sensing capabilities. 
In Fig. 6(b) we plot Zopt versus SNRc with p as a parameter. 
For a given SNRc, we see that it is more profitable to increase 
the size of the cluster when the process is highly correlated 
than to increase SNRpfs-mac(0- 




(b) p as a parameter. 


Fig. 6. /opt versus SNRc with e = 2, /3 = 1 and n = 100. 


The optimum cluster size ^opt allows to redefine the idea 
of spatial coherence regions (SCRs). Typically, a SCR is 
an inherent property of the process to be sensed, uniquely 
determined by the spatial correlation between its samples |j^. 
In this work, guided by the error probability as an indicator of 
the ultimate performance, we obtain, that a SCR is not only 
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defined by the spatial coirelation of the process, but also by 
critical parameters as SNRm, transmitted energy E, path loss 
exponent e, and bandwidth /3 available at each cluster. 


B. Random Network 


In this section, we compute the performance of the statistics 
when the sensors of the network are spatially distributed as a 
homogeneous Poisson Point Process of intensity A > 0. 
In this case, the total number of sensors and the number of 
sensors per cluster follow a Poisson distribution with mean 
XuAq and XIAq, respectively. For simplicity, we continue 
using the assumption that the clusters are squares of area 
IAq and the region to be sensed has a total area uAq. We 
set A = 1 and Aq = 1. Therefore, the average number of 
sensors per cluster coincides with 1. We also assume that the 
spatial process is isotropic meaning that the ACF depends 
only on the distance d between the points for which it is 
evaluated. Some examples of possible isotropic ACFs are the 
equicorrelated process presented previously, the exponential 
ACF Rs{d) = or the hyperbolic ACF Rs{d) = 

The exponential ACF was found to be a good model for the 
shadowing effect in a urban cellular environment p8] and it 
will be used here. 

The precoding strategy used for an arbitrary spatial distri¬ 
bution of the nodes is AFS-PAC whose gains are defined in 
(|^. For this strategy, the signal-to-noise ratio for the fc-th node 
at the i-th CH is defined as SNRAFs-PAc(fc) = 

Replacing (|^ in it gives 

cvTD SNRcSNRm 

SNRafs-pac(^) = 


l+alla 


SNRc + (1 + SNRM)||atfc — 


We compute the average miss error probability as a metric 
of performance, defined by the average over the PPP $, 
Pm = E(Pm($)) subject to Pfa($) = a for each realization 
of the PPP. This average is computed using the method of 
Monte Carlo and the theoretical expressions found in Corollary 
[2 (labeled LCS, FCS and ED in Figs. and |8]l. On the 
other hand, we also generate the sensed process and the 
corresponding noises using again the method of Monte Carlo 
to validate the theoretical expressions derived in Corollary [T] 
(labeled LCS-MC, FCS-MC and ED-MC in Eigs. |7] and [8]l. 

In Eig. |7] we plot the miss probability as a function of 
the average number of nodes per cluster I with the param¬ 
eters indicated in its caption. Notice that the results for a 
random network with exponential ACE are similar to the 
equicorrelated process case. Considering the LCS strategy, in 
spite of decreasing the effective signal-to-noise ratio at each 
CH SNRafs-pac(^) for most of the sensors when I increases 
up to 10, the miss error probability decreases because of 
the intra-cluster correlation is better exploited. If the cluster 
size increases further, the path loss effect inside each cluster 
dominates and the performance is deteriorated. The LCS 
strategy is close to the bound provided by the ECS scheme 
and the gain compared with the ED is considerable. 

In Eig. we plot the miss error probability as a function 
of SNRc with the parameters indicated in its caption. Notice 



Fig. 7. Pm versus I for an exponential ACF with p = 10, n = 144, 
= 10-2, SNRm = 0 dB and SNRc = 12 dB. 


again that the results are similar to the closed-form expressions 
obtained for the equicorrelated process. The LCS strategy is 
close to the bound provided by the ECS scheme and the 
gain respect to the ED strategy is considerable for the entire 
shown range of SNRc. These Monte Carlo curves also validate 
the results of Th. [T] where the central limit theorem was 
applied to compute the error probabilities, obtaining good 
approximations even for a low number of sensors n. 



Fig. 8. P" vs SNRc for an exponential ACF with p = 10, n = 144, 
I = 18, m = 8, Pfa = 10—2, and SNRm = 0 dB. 


VII. Conclusions 

In this paper, we developed a model and appropriate tools 
for evaluating the performance of a distributed detection 
scheme that makes a trade-off between transmission energy 
and exploitation of the spatial correlation of the process. 
We characterized this trade-off and found two regimes of 
operation: an energy limited regime, and a correlation-limited 
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regime. We showed that the proposed LCS strategy allows to 
obtain excellent performances, which are close to the ones 
obtained with the use the global network-wide correlation 
information. The use of the spatial correlation is of great 
importance in WSNs where the captured measurements have 
low signal-to-noise ratio and also when the process to be 
sensed presents moderate to high correlation. 


Appendix I 

Error Probability Approximations 

In this appendix we derive an approximation for both the 
false alarm probability and the miss-detection probability 
following a line reasoning similar to p9) and using a 
version of the central limit theorem (CLT) for a sum of 
independent but not identically distributed random variables 
given by the following lemma. 

Lemma 2 (Essen CLT). Let Xk be mutually independent 
random variables such that 


E(a:fc)=0, E{xl)=al, E{\xk\^) = pk- 


Let = (Tj -I-h cr^, r„ = Pi -I-h p„ and denote by 

the cumulate distribution function (CDL) of • •-|-x„)/s„. 

Then, for all u and n, 

\Fn{u) - $(u)| < 

where $ is the standard normal CDL. 


Proof: See 130 p. 544, Th. 2] 


A. Proof of Theorem 

We begin with the computation of P(T„ > r„) when > 
E(T„). Let Qk be the PDL of yk,s, the tilted version of yk 
with PDL pk, such that qkivk.s) = . It is 

straightforward to check that the first and second derivatives 
of p.k{s) satisfy /ifc(s) = E(yfe,s) and flk{s) = YaT{yk,s), 
respectively. Define = yi,s + ■ * * + yn,s- The PDF of ^ 
is 


P n 

Qx^fs) — / 11 Qk(,yk,s)t^yk,s 
k=l 
. n 

= / l[pk{yk,s)e^^^’^-'^>‘^^Uyk,s, 

k=l 

where the integration domain TZ is dehned as TZ = 
{(t/i, G K” : 2 / 1 ,s H-h yn,s = Therefore, 


Qrits) 


yj f[pk{yk,s)dyk,sje''*^ 


where Pt is the PDL of T„ and /r^(s) = J2k=i Pkis) is its 
LMGL Thus, r„ ^ is the tilted version of T„. Then, 

/ OO 

PT{f)dt 

/ OO 

(30) 

poo 

J Lillis) 


(31) 


The first factor in (30 1 is equal to the Chemoff bound if s 


is constrained to be positive. Because the LMGL /ry(s) is 
a convex function, fris) is an nondecreasing function. By 
assumption, > E(T„). Thus, the value of s that solves 
Tn = At(®) satishes s > 0 and we can use this equality to 
replace r„ in ( [3T] i. 

The Chernoff error exponent dehned as 
lim„_>oo ^(sAt'('®) ~ Pri^)) asymptotically tight, i.e., it 
is equal to the error exponent lim„_>oo — ^ logP(T)i > r„). 
Therefore, we need to approximate the integral which is a 
sub-exponential factor. We let Ft be the CDL of Qt and 
integrate by parts to obtain 

poo 

poo 

d Pt(s) 

y/I?^^“(GT(u)-Gr(0))e-«V^“du, (32) 


= 5 


where in the second equality we have defined the following 
zero-mean unit-variance random variable: Un = _ 

jygjj P Applying Lem. 2 with 

V4t('’) I 

Xk — yk,s — Pk{s), it can be shown that s\ = jl'ffs) = 
0{n) and r„ = 0{n). Therefore, Gt{u) — Gt(0) — ^{u) — 
<h(0) + 0{^), where 4) is the standard Gaussian CDL. The 
exponential factor in ( [32) i suggests that for sufficiently large 
n we only need a good approximation of ^(m) — 4?(0) around 
0, and only within a small fraction of the standard deviation 

of Tn, J'Pt{s). Thus, it is sufficient to use the following 
2 

bounds: 1 — \ < e~^ <1 which produces, if u > 0, 
^ i ^ ^ Then, it easy to show 

that 

v'27rs2/ij(s)) Vv^ 

Replacing this in ( [31] ) produces ( |2^ , the hrst result of the 
theorem. The computation of P(T„ < r„) is similar and it is 
omitted. 


B. Proof of Corollary [^ 

In order to apply Thm. [T] the Gaussian quadratic form ( |25| ) 
can be expressed as the sum of n' independent variables. Then, 
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we only need to evaluate its LMGF fii{s) = logE^ 




its first and second derivatives are: 


— log det ■} 

(33) 

Mo(so) 

diis) = tr{[(5,_„-P„-)"^ - s4/]“^}, 

(34) 

Po{so) 

Pi(s) = tr{[(5j,„/P„-)“^ - s/„/]“^}. 

(35) 

Replacing these expressions in (231 and ([24| completes 

the 

Po{so) 


proof. 


Mo(so) = - jlog (l - So 


r(i+(i-i)p) 
i+r(n-(z-i)p) 

- n (/3 - ]-) log (l - so i+[^\7-p) 
r(i+(i-i)p) 




iP-\) 


r(i-p) 

i+r(i-p)-sor(i-p)> 

(r(i+(;-i)p))^ _ 

:i-i)p)-sor(i+ii-i)p)y 

(r(i-p))" 

(i+r(i-p)-sor(i-p))=- 


(36) 


(37) 


(38) 


Appendix 11 

Performance for the equicorrelated process: 
PROOF OF Lemma[T] 


Now, the mean of the statistic Tlcs under "Hi is mi = 
El (Ties) = nr{p + /3(1 - p)). Then, ^o(sS) = r„ = mi 
implies that Sq = 1. Replacing it in ([36]l-([38| and using Th. 
we obtain ( |27| ). 


A. Pfa for the LCS strategy 


We use Th. 1 to compute Pf 


fa,LCS- 


We hrst evaluate the 


LMGF of the statistic Tlcs iu ( [20l i under "Hq und its hrst and 
second derivatives given by (p3])-([T5|) in Corollary [T] where 
P / — r / 

J- n.' — -L r).' • 


The covariance matrix of the equicorrelated process is 
a circulant matrix and, therefore, it is diagonalized by the 
I X I DFT matrix Fu and its eigenvalues are the DFT of 
its hrst column, i.e., = FuDgiF^ where D^i = 

diag(l + (Z — l)p, 1 — p,..., 1 — p). Considering the PSD 
of the equicorrelated process, the frequency vectors used in 
the precoding matrix of the PFS-MAC scheme given in Def. 
1^ have the indexes jk> = k' — 1, k' = 1,..., T. Then, the 
covariance matrix in is 


Si,P = 7 o-Ds./' + {-Jocrl + ai)Ip, 


where Dgj' is a T x I' diagonal matrix, = diag(l + (^ — 
l)p, 1 — p,..., 1 — p). Notice that Si_// is a diagonal matrix 
because the covariance of the process is a circulant matrix 
and the precoding scheme is based on the DFT matrix. The 
covariance matrices of 2 ; under T-Lq and Tfi are given by ( [T5| ) 
and they result 

^0,n' = ilo^v 3“ 

Si,n/ = 72 cr 2 (l_p)/^, + Fif,)li{Im Fip) 

+ (lo^v + 

= 7o'^s((l ~ P)In’ + P'i-m C) Ep) + (7oCr^ + 


where 0 is the Kronecker product and the I' x I' matrix Ep 
is dehned as Ep = diag(/, 0,..., 0). Then, 


^0,n' 


V , — T 

X n' — Xr7 


’ diag 


Ipdl 


Igdp 


ll{d^+al)+al ’ ■ ■ ■ ’ ^l{dp+al)+al 


where di = (jf(l + {I — l)p) and dk = crlil — p) if k = 
2 ,... ,1'. Using ((3^-((3^ of Corollary [T] and denoting F = 
SNRpps-MAc(0’ '^he LMGF of the statistic (p0|) under Hq and 


B. Pfa for the FCS strategy 

The statistic Tpcs is an LLR, for which the LMGFs under 
"Ho und Tdi satisfy 


Po(s) = /il(s - 1), 


(39) 


where Pi{s) = Ei(e®^‘'“). The direct computation of po('S) 
implies to obtain firstly the inverse of Si_„/, which is not 
straightforward. An easier way is to compute pi{s) for which 
we only need to invert the diagonal matrix So^ri' and then, to 
use (39 1 to evaluate po{s). The expression for pi(s) and its 
first and second derivatives are: 


pi(si) = 

- Si log (1 + r(l + (n - l)p)) - Blind - 1) log (1 + r(l - p)) 

- log (1 - Sir(l + (n - l)p)) - (n/3 - 1) log (1 - sir(l - p)), 
pi(si) = 

- log (1 + r(l + (n - l)p)) - ind - 1) log (1 + r(l - p)) 

r(i + (n-i)p) , , ^ . r(i-p) 


l-sir(l + (n-l)p) 
Pl(si) = 

(r(l + (n-l)p))^ 

(1 - Sir(l + (n - l)p))2 


+ {nd - 1) 


+ {nd - 1) 


i-sir(i-p)’ 

(r(i-p))^ 

(i-sir(i-p))2' 


Tpcs = Ei(rpcs) = /ii(si) = Po(so) implies that si = 0 and 
So = 1- Using this in ( |39] l and replacing all in ( [2^ , we obtain 
( |28] ). Notice that this result coincides with the one in ^n\ if I 
is replaced by n. The difference between both FCS and LCS 
strategies is that the first one contemplates the full correlation 
of the process. 


C. Pfa for the ED strategy 

The energy detector Tpo can be normalized by 7 o<t^ + cr^ 
without modifying the performance: Ted = .^ 2 ^™ ^2 ■ Its 
LMGF under "Ho and first and second derivatives are: 


Mo(so) 

Ao('So) 

Po{so) 


-nlog(l - So), 
n 


l-so 

n 

(1^^' 


(40) 

(41) 

(42) 
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Ted = Ei(Ted) = n{l + F) = /to(so) implies that sq = 
Using this in ([40)i-([42)i and replacing all in ( |2^ results in ( |29] l. 
Notice that the same result can be obtained by replacing I by 
n in (271 and taking p = 0, given that the ED strategy does 
not consider the correlation of the process. 
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